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INTRODUCTION 

Under Grant NSG- 1-206, members of the Faculty in the Department of 
Aeronautics and Astronautics at MIT, together with selected 
undergraduate students, have undertaken four projects as part of the 
Undergraduate Research Opportunities Program at MIT. The research under 
this grant is coordinated with specific Langley personnel, but not among 
the individual investigators. The remaining part of this report will 
consist of individual summaries by each of the four investigators. For 
convenience, the investigators are listed in alphabetical order. 
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PROGRESS REPORT: HOT WIRE IN LOW REYNOLDS NUMBER FLOW 

Mark A. Kolb 
Prof. E. E. Covert 

Project Objective 

The intent of this project is the study of the low Reynolds number 
flow phenomenon of periodic vortex shedding in the wake behind a 
cylinder, applying the hot wire anemometer technique of flow 
measurement. This work is concerned primarily with the downstream 
diffusion of these shed vortices. An evaluation of the performance of 
the hot wire at low Reynolds number will also be considered. 

Progress to Date 

This work began slowly last spring with the first faltering steps 
at becoming familiar with the facilities to be used at the MIT Gas 
Turbine Laboratory, where the work is being done. At this time, the 
concepts behind hot wire measurements, as presented in the literature, 
were also investigated. 

Work in the summer began with the correlation of two sets of 
pressure measurements. The vortices are produced in a free jet 
(one-inch diameter) which flows from a tank connected to an oil-free 
compressor. Pitot tube measurements of the jet flow velocity were 
correlated with pressure tap measurements from the wall of the tank 
from which the jet flows. This enabled the hot wire calibration to be 
performed with respect to the pressure tap measurements, instead of 
against pitot tube measurements; simultaneous pitot tube and hot wire 
measurements might have been mutually inaccurate due to interference 
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effects inside the small jet. 

Progress quickened mid-summer when the opportunity arose to begin 
constructing hot wire probes, and the work of constructing, 
calibrating, and repairing probes occupied most of the remainder of the 
summer. A sketch of the probe design is presented in Figure 1 . Gold 
electrodes are soldered to a pair of broaches (round, dimensioned 

"0.0003" x 0.013""). The broaches are then inserted through a short 

' 

length of ceramic tubing. The ends of the broaches, which protrude 
about half an inch from the tubing, are then bent in half to a 90° 
angle, 2.5 :m-diameter platinum wire, Wollaston-drawn, with silver 
coating, is used for the "hot wire" itself. The silver coating is 
etched away from a few millimeters of wire, using nitric acid. This 
etched length is then soldered across the ends of the two broaches. 
Present Work 

Work has continued this fall using three of these hot wire 
probes. A DISA 55D01 Anemometer is employed to process the probe 
signals. A sketch of the present apparatus is shown in Figure 2. 

The calibrated hot wires are being used to determine optimum 
conditions under which the wake vortices may be observed. Conditions 
such as flow speed and cylinder diameter have been varied while 
observing the output of the hot wire on an oscilloscope. 

At first, several problems were encountered. Not considering the 
effects of the vortex's passage, sinusoidal output was originally 
sought. Figure 3, showing a typical wake velocity-time profile, 
indicates that, indeed, the flow is periodic - as expected - but it is 
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not sinusoidal. The passage of a vortex is indicated by the sudden 
bursts in output voltage. It appears that there may be an underlying 
sinusoidal signal (it is difficult to discern whether or not this is 
actually just an illusion created on the rapidly changing oscilloscope 
screen), but it is not likely that this is due to the vortex shedding. 
On the other hand, it does not appear in the trace of the free jet hot 
wire measurements (i.e., without the cylinder). 

Originally, cylinders of too-small diameter were used, but the 
vortices were too difficult to locate. Present measurements (Figure 3 
included) are for a 4.2 mm cylinder in the one-inch jet. Location of 
the cylinder in the jet (i.e., outside the zone of potential flow) was 
another hindrance to good results. A third problem concerns the flow 
velocity. To improve the vortex formation, flow speed was increased. 
The 0.96 V average output voltage in Figure 3 corresponds to a flow 
velocity on the order of 1 .4 m/sec, for a Reynolds number based on 
cylinder diameter of about 380. This is well within the range for 
periodic vortex formation. If, however, the upstream speed is 
considered, a Reynolds number of almost 1800 is observed, which 
suggests a turbulent wake. Perhaps, then, the results presented in 
Figure 3 are not optimum. 

Intentions 

Upon establishing a working ability to distinguish the wake 
vortices, it is next desired to begin quantifying the wake parameters. 
The first measurements to be taken will concern the dimensions of the 
wake - in what region may the vortices be found? How far above and 
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below the cylinder? How far downstream? This can be accomplished 
simply by moving the hot wire probe and measuring the displacement with 
rulers. Preliminary results suggest that, at least for higher flow 
velocities, the edges of the wake are quite distinct. 

The next quantity to be determined will be the speed of diffusion 
of the vortices downstream. Measurements of the frequency of the 
periodic hot wire signal at different locations in the wake will 
describe the velocity profile of the shed vortices. Another quantity 
to be measured might be the strength of the vortices, but a method for 


doing so has yet to be considered. 
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connection to 
probe holder) 


Figure 1. Sketch of hot wire probe design. 
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Figure 2. 


Sketch of present experimental set-up. 
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Figure 3* Sketch of sample oscilloscope 

trace of the anemometer output. 





Flow Behavior in the Wright Brothers Facility 


Investigators: Prof. Judson R. Baron 

Mr. Stephen Genn 

NASA Technical Monitor: Dr. Robert Kilgore 

At various times during the last several years test results 
obtained in. the Wright Brothers 7.5' x 10' wind tunnel have led to 
speculations as to the flowfield uniformity. Such test considerations 
are typical and always proper, and in combination with numerous 
alterations that have been made recently to the tunnel circuit did 
suggest that the flowfield should be reexamined. 

The objective was seen to be twofold: 1 ) A brief examination of 

back sections of the circuit to establish whether or not gross flow 
deviations were present at corners, turning vane regions, etc., and 2) 

A calibration of the test section itself. 

Both of these tasks were undertaken during the year and both the 
procedures and results of the project are described in the attached 
report (Flow Behavior in the Wright Brothers Facility, Stephen Genn, 
WBWT-TR-1 187, September 1983). In essence rather substantial 
calibrations were obtained with total and static pressure 
instrumentation over large portions of two test section vertical 
planes, and five speeds in the range 35 to 125 mph (50 - 185 fps). The 
useful central core flow proves to extend over approximately 0.75 of 
horizontal and vertical spans of the elliptical test section, with 
dynamic pressure variations over the entire section being within 1%. 

The report is attached as Appendix A. 
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Final Report on Development of Software for 
Digital Signal Processing 

Investigators: Professor W.L. Harris, V. Del Gaudio, 

R. Edmonds, and S. Zachary 

NASA Technical Monitor: Mr. D.R. Hoad, LSAD 

The research efforts of this project were concentrated on the 
development of auto- and cross-covariance and of auto- and cross- 
correlation software. This Fortran code was written from a time and 
frequency domain analysis. A second development consisted of a code to 
compute the power spectrum. Both codes have been debugged and are 
available on IBM compatible 8" flexible disc. 


- -r*#*** 
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Title of Task: Numerical Studies of Some Aspects of the 0.3 Meter 

Cryogenic Self-Streamlining Wind Tunnel 
Principal Investigator: Professor Earll M. Murman 

Student: Jonathan Webb, Class of 1983 

Final Summary of Progress : 

This task involved developing a numerical simulation code to 
assess residual errors on airfoil loads due to imperfect shaping of the 
self-streamlining wall for the 0.3 meter cryogenic wind tunnel. Causes 
of imperfect shaping could result from pressure loadings on the wall 
and the finite number of jacks available to shape a wall representing a 
freestream streamline. It should be noted that previous development 
work for the wall included an estimate of such errors. However, the 
estimates had been made without using transonic computational methods. 
The current task involved writing a simple numerical model for the wall 
structure, modifying an existing transonic small disturbance code for 
the free air and self streamlining wall conditions, coupling the two 
codes together, and performing numerical simulations to quantify the 
residual errors. 

A detailed summary on the progress of the research is contained in 
the attached report entitled "Adaptation of Two-Dimensional Transonic 
Analysis Code, TSFOIL, for Use in Modeling Adaptive Wall Test Sections" 
authored by the student (Appendix C). The structural code was written, 
the TSFOIL code modified, and the two codes hooked together. 
Unfortunately, by the time the student completed his undergraduate 
studies in May 1983, the calculations were not converging. It is not 
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known if there are errors in the formulation or the coding. This proved 
to be a rather challenging task to carry out within an undergraduate 


research project environment. 



Title of Task: Induced Drag Predictions for Aircraft Configurations 

Principal Investigator: Professor Earll M. Murman 

Student: Richard Shapiro, Class of 1983 

NASA Technical Monitor: Edward Waggoner 

Final Summary of Progress : 

This task was terminated shortly after it was initiated due to the 
unavailability of the student. No funds were expended on this task. 
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Title of Task: Calculating C-Grids with Pine and Embedded Mesh Regions 

Principal Investigator: Professor Earll M. Murman 

Student: Bernard Loyd, Class of 1983 

NASA Technical Monitor: None 

Final Summary of Progress : 

This task was initiated during the summer of 1983 due to the 
availability of the student on rather short notice and the timeliness 
of the research topic. Research by a PhD student, Mr. William Usab, 
demonstrated the attractiveness of using embedded grids for airfoil 
calculations modeled by the Euler equations. Mr. Usab's work had been 
done exclusively using O-grids. It was desirable to extend these 
calculations to C-type grids both for Euler, and eventually, the 
Navier-Stokes calculations (see Appendix C). 

During the summer of 1983, Mr. Loyd developed a C-grid generation 
code for an arbitrary airfoil shape. The method followed the 
Jameson-Caughey approach of a parabolic mapping followed by shearing 
and stretching transformations. A working and documented code was 
developed during the summer on a VAX 11/750, and in the fall Mr. Usab 
utilized the code for Euler and preliminary Navier-Stokes calculations. 
Comparison of the Euler results with the previous 0-grid calculations 
demonstrated significant improvements in rates of convergence and 
required levels of artificial viscosity. 

During the fall of 1983, Mr. Loyd extended the C-grid code to 
include the capability of embedding any number of refined grids within 


the solution domain. A detailed description of the work is included in 
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the attached report entitled "Calculating C-Grids with Fine and 
Embedded Mesh Regions" authored by the student. Mr. Loyd plans to 
continue working in this general area for his Masters thesis research. 
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Closing Remarks 

In the last two years, eight undergraduates in the Department of 
Aeronautics and Astronautics were provided with opportunities to study 
problems of importance to several research tasks. The problems varied 
widely, from writing programs for computer aided experimentation to 
attempting to measure the decay rate of the vortices that make up the 
Karman vortex street. This latter issue is related to interpretation of 
certain flow visualization processes, i.e., does one have a long vortex 
street or just the marking points where the decayed vortex was? 

In any case, the faculty feels the educational opportunity offered 
by this grant was great and that these students profitted greatly from 
their experience. As an experiment in a cooperative educational 
venture, we feel it was successful. 
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1. Introduction 


It has become increasingly apparent that a reexamination of the flow 
characteristics in the low speed Wright Brothers Facility (WBF) is of some 
importance in view of recent improvements in the precision of the data 
acquisition system, in particular, the existence of local regions of 
separation, if any, in' back portions of the circuit, and possible related 
unsteadiness, are of interest. 

During April 1983 a preliminary qualitative experiment was performed 
to examine the nature of any problem. The test was performed in the up- 
stream cross leg preceding the stilling section (marked Section D on 
Figure 1-1) . Observations from that initial experiment did indicate some 
unsteady air flow problems in the cross leg, and thereafter the test region 
(Section A) was calibrated quantitatively. The intent was to learn some- 
thing about the effect of upstream intermittent behavior flow on the test 
section flow, as well as to provide an extensive calibration as a standard 
for the effects induced by future alteration of the tunnel. 

Distributions of total pressure coefficients were measured first at 
one cross-section plane of the test section, namely the model station 
(Fig. 1-1) . Data were obtained for several tunnel speeds. The reduced 
data yielded an unexpected distribution involving larger pressures along 
the inside wall. It was conjectured, for example, that a rotation of the 
entire channel flow may exist. Consequently, a second calibration was 
performed in a plane located seven feet downstream of the model station. 

The following sections describe the test conditions, procedures and 
results for the present unmodified tunnel circuit. 

2. Preliminary Observations, Back Leg, Section D 

On April 8, 1983, the initial observations were made in the Wright 
Brothers Facility in Section D (Fig. 1-1) . Although a grid pattern 
of tufts was conceived to aid in cataloguing observations, only the 
flow in the very central plane of Section D was observed (Fig. II-l) . 

The observed plane, however, was divided into a grid system as shown 
in Fig. II-2 . 
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Fig. II-l. Plan view of Section D. 


Fig. II-2. Cross-sectional 
view of grid at plane 
indicated in Fig. II-l. 


A ten foot conduit pole with six two-foot yarn tufts spaced at six-inch 


intervals was used to observe the flow (Fig. II-3) . All grid areas were 
within reach excepting #2 and #3 (Fig. II-2) . The tunnel was run at a 
constant 40 mph (in test section) and photographs were taken to record 
phenomena of special interest. 
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Fig. II- 3. Photograph 
of tuft pole in 
Section D, looking 
towards inside wall. 



OBSERVED: 
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Center region grid sections #6, 7, 10, 11: ! 

All flow in these sections was observed to be steady. 

Boundary grid sections #4, 8, 12, 16, 15, 14: ' j 

All flow, exclusive of a four to six inch boundary layer, was 
steady. [ 

Inside wall grid sections 1, 5, 9, 13: 

The flow did not appear to be turbulent in this region. However, 

J 

there was a surging of sorts corresponding to., a cyclical growth and decay 


of the flow velocity. The oscillatory velocity had no simple time de- 
pendence. The cycles were periods of the order of five to ten seconds 


(Figs. II-4 and II-5) . 

Fig. II-4. Double exposure 
showing the decaying flow 
phenomenon on the inside 
wall. Photo taken from 
grid section 11; view is of 
tufts on the border between 
sections 5 and 9. 




, Fig. II-5. Tufts in ri 

normal steady flow of 1 f 

air. Photo taken five 
seconds after exposures 
shown in Fig. II-4. j 

n 


< 
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3. Calibration of Test Section, Model Station 

Total pressures were measured at the model station plane of the test 
section. Some static pressure measurements were also taken. 

The tunnel was operated at five nominal velocities: 35, 50, 70, 100, 
and 125 mph. A device resembling a long rake with -21 pitot probes at 
four-inch intervals was moved through nine positions in one plane of the 
test section at the model station (Figs. III-l, III-5) . Consequently, 
during this calibration there were 9 x 5 -45 test runs. 



In addition, the following were measured: total pressure at the up- 
stream entrance to the test section, the average of four static pressures 
located at flush wall taps on and spaced evenly in a circumferential ring 
around the test section, and both total and static pressures at three 
pitot-static probes mounted in the test plane (Figs. III-2a and III-2b) . 
The pressure at each such tap was measured during each test run. 

The rake was held in place on two steel unistrut P-400 bars with 
C-clamps (Fig. III-3) . 
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Pitot probe j-fsr section 



Fig. III-2a. Side view of 
test section showing 
upstream reference pitot 
probe and the calibration 
planes . 







STATIC TAPS : Four taps around i 

section as indicated; their 
average pressure provided the 
single reference static pressure.^ j 

PITOT-STATIC PROBES: Three 
probes in test plane. 


Fig. III-2b, Cross-section of the model station test plane, showing the 
positions of pitot-static probes and reference static taps. 


The measurement system consisted of a Type D, solenoid operated 48-tap 
scanni valve in combination with a model 237 ±0.5 psi Setra pressure trans- 
ducer. Plastic tubing was used between the scanni valve and the taps. A 
Baratron unit was also used to measure the upstream pitot pressure. Both 
Baratron and Setra units measured pressure differences relative to the 
reference static pressure (Fig. III-2b) . The Setra unit was mounted on the 
back of the rake, with remote input and all output in wire and tubing 
entering through a hole in the floor (Fig. III-4) . 



Fig. Ii; 
rake poi 
tunnel. 


Fig. Ii: 
scannivi 
transdui 
rake. 


run. Fol 
d at each 








1 Reference total pressure upstream of entrance to test section i 

(Fig. III-2a) 

2 - 22 Total pressure at 21 pitot tubes on movable rake (Fig. Ill- 3) 1 


24, 23 Static and total pressures at pitot-static probes (Fig. III-2b) 

in measurement plane, near top, center. 

25, 26 Static and total pressures at pitot-static probes (Fig. III-2b) j 

in measurement plane, towards outer wall, i.e. on left facing 
upstream. 

27, 28 Static and total pressures at pitot-static probes (Fig. III-2b) ' ; 

in measurement plane, towards inner wall, i.e. on right facing 
upstream. 


The scannivalve position was controlled and the data recorded with 
the aid of a computer program (WBWT1, prepared by F. H. Durgin and 
J. Biller) which generated as output a pressure coefficient for each tap. 

The total pressure (relative to reference static) measured at tap 1 was 
used to normalize all pressure measurements. The data were corrected for 
instrument offset and drift, and plotted as functions of position. Finally, 
contour maps were prepared to show loci of constant normalized total pressure 
parameters in such test planes. A map was drawn for each velocity. For a 
more detailed explanation of the numerical data reduction, see Appendix I. 


Computer Output and Data: 


r*\ 





n 

I 


During the calibration test a distinct file was generated as output 
for each test run. Since there were nine rake positions and five speeds, 
a total of 45 output files were generated. Each file for this first data 
plane (at the model station) is named TUNCAL.abc, including a three-digit 
descriptor abc. The first digit, a, of the descriptor indicates the 
nominal wind speed: 

1 35 mph 

2 50 mph 

3 70 mph 

4 100 mph 

5 125 mph 

The second digit, b, of the descriptor provides for specifying future 

calibrations. In this case, the second digit is always 0. The third 
digit of the descriptor, c, designates the position of the rake in the 
test plane (Fig. III-5; photographs appear in Appendix II). 
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4. Calibration Data; TUNCAL. 101-509 
Appendix III, Part A, contains data from the calibration carried out 
in the model station plane. The first three pages in the Appendix are 
tables of computer output. TheCp n entries in each column are a series of 
normalized pressures; each column of computer output corresponds to a 
single test run. A more detailed explanation of the computer output and 
the subsequent numerical corrections for instrument shift, etc. , are 
described in Appendix I. 

Following the reduced data are six pages of graphs (App. Ill, Part B) 
showing normalized total pressure parameters as a function of position in 
the test plane. Positions 1 and 2, positions 3 and 4, and positions 6 and 
7 share pages respectively (see Fig. III-5) . There are five graphs included 
on each page, one for each velocity. 

Finally, there are five contour maps, one for each speed, showing 
contours of constant total pressure parameter (App. Ill, Part C) . These 
maps are based on interpolations from the graphs of total pressure parameters 
with respect to position. 

5. Calibration of Test Section, Second Measurement Plane 

This series of test runs was identical in purpose to the series at the 
model station (Section 3) excepting that the rake and the pitot-static probes 
were positioned in a second test plane located seven feet downstream of the 


^' V; '£T 
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model station (Fig. Ill-2a) . Test runs were performed at only three speeds, 
and, unfortunately, the fourth and ninth rake positions were omitted due to 
a tunnel accident. 

The output data file names for this series of test runs were of the 
format TUNCA2.abc including a three-digit descriptor following. As before, 
the first digit of the extension indicates velocity and the last refers to 
rake position (Fig. III-5) . For this series, the second digit is always 4. 

Note that the upstream total and the four wall static reference taps 
(i.e. taps 1 and the average on 48) were not moved for the second test series. 
Those measurements refer as before to the same tunnel locations 
(Fig. III-2a,b) . 

The data for this second test series appears in Appendix IV, Parts 
A, B and C. 

6. Conclusion 

Parts B and C of Appendices III and IV indicate the following: 

A. Dynamic pressure distributions across test section are in most 
cases within ±1% over half the span, and for speeds in the range 50 to 100 mph. 

B. A bias towards higher/lower speeds on the inner/outer sides (i.e. 

right/left in figures) is apparent near the horizontal midplane (position 

pairs 1, 2; 3, 4; and 6, 7) . The contours drawn in Part C illustrate the 

»» 

relatively accelerated regions (larger Cp n ) more clearly. 

C. The vertical survey (position 9) shows quite flat variations 
along the vertical midplane. In combination with the horizontal surveys, 
a slight antisymmetry about the vertical midplane is present but with a 
useful central core flow representing about 0.75 of the horizontal/vertical 
spans . 
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APPENDIX I 
DATA REDUCTION 


Two measurement systems were used during the program. A Setra pressure 
transducer connected to a scannivalve was used for all taps. A 
Baratron unit also was used for the upstream pitot (tap 1) and its indicated 
pressures were used as a basis for normalizing all other readings. Both 
Setra and Baratron units were referenced to the average static pressure in 
the model plane (first calibration) as determined from the four wall taps 
around the tunnel boundary (recorded as tap 48) . 

Define: 


p. total pressure at upstream pitot (tap 1) 

t REr 

p static pressure average based on four wall taps in 

model plane (tap 48) 

p pressure, total or static, at tap n 

(see description, p 8) 

q REF reference dynamic pressure (p fc ^ - P g ^ ) 


The pressure transducers (Setra and Baratron) are referenced to indicate 

increments above p . Thus 

s REF 


P BARA P t REF P s REF ^ q REF 
(P SETR } P n ~ P s REF 


(for tap '1) 
(for tap n) 


n 

if the instrument systems were perfect. Actually there are errors 

due to initial offset and zero drift during the measurement procedures. 

Say 


( 2 ) 


(p ) 
^s REF' 


= P 


Apparent 


s REF 


+ Ap + Ap 

^OFFSET *DRIFT 


(3) 


in which the offset and drift contributions may differ for different trans- 
ducers. Then the recorded values would be 


P = Q 

BARA ^REF 


^OFFSET + AP DRIFT^ 


^■ P SETR^ 


n 


SETR' 


BARA 

'OFFSET + Ap DRIFT-^ 


SETR 


( 4 ) 


The value of (P SETR > pri°r to a data scan establishes < A P 0FFSET ) SETR since 
the initial drift is zero by definition. I.e. from eq.(4b)* 


^ P SETR.^ 481 " ^ Ap 


OFFSET' SETR 


( 5 ) 


*This initial .zero reading is not included in the tables. 
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and has been used to remove that tare from all of the Setra data. [Subscript 
481 refers to the "initial" tap 48 reading.] The drift contribution to the 
Baratron readings is known to be small and it is assumed here that 


^ P DRIFT^BARA ° 

Now eq. (4) , with (5) and (6), becomes 

P BARA ~ q REF “ ^ P OFFSET^ BARA . [= P< BARA 
^ P sf.tr) - ^ P SETR^ ~ ^ Ap nRTF r r)qF'rp (“^P 


■ DRIFT'SETR 


SETR^n (P SETR } 48I^ 


n n 

and we define an uncorrected, normalized pressure coefficient to be 

(p" ) 

^ SETR'n 


Cp n = 


fp" ) 

^ BARA' 


( 6 ) 


(7) 


(8) 


Cp n is tabulated in Parts A of Appendices III and IV for the data that was 
recorded for the two calibration planes. 




SETRA Drift Correction 

The last read tap 48 outputs shown in the Appendices indicate the total 
drifts during a given scan, i.e. from eq. (7) 

4Pt ~ ^ Ap nRTFT^ ! 


DRIFT-' SETR, TOTAL 


(9) 


The correction will assume that drift is linear with time. Since the scan 
is also sequentially linear over all N(= 29) taps. 


^ AP DRIFT^SETR,n N ^~ P 


Then 


Cp = 
n 


^ P SETr) n 


BARA 
n 


= Cp + 
n 


SETR 48^ 
^ P iirtft) 


DRIFT'' SETR 


BARA 


= Cp n ~ i Cp 48 
and is a value corrected for drift. 


( 10 ) 


( 11 ) 


-n 
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Baratron Offset Correction 


From Eqs. (11), (7) and (2), for tap Is 
^- p setr) 1 1 


Cp. = - 

i q. 


REF 


OFFSET^ BARA 




lAp, 


7 


( 12 ) 


OFFSET^ BARA 


1 REF 


and the offset correction is implied by the difference from unity of Cp^ 
Eq. (11) may now be written as 


CP 


' ( p SETR^ n o' 


n 


CP, 


(13) 


REF 


or 


fp 1 Cp' Cp - - Cp 
" _ u SETR' n _ n _ n N . 

Cp 1 


% = 


48 


(14) 


REF 


1 C %- 5 <*48 


is a fully corrected. pressure . coefficient corresponding to a ratio of the 
"perfect" differences defined in Eq. (2) . The normalized, corrected, 
pressure coefficient Cp^ has been found from the tables in Appendices III 
and IV and is plotted and contoured (after interpolation) in Parts B and 
C of those Appendices. 
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APPENDIX III, Part A 

Uncorrected Calibration Data for TUNCAL.abc Files, First (Model Station) Plane 


FILE. abc 
Vel. ,mph 

.101 

37 

.102 

34.7 

.103 

34.5 

.104 

34.5 

.105 

34.5 

.106 

32.8 

.107 

33.3 

• 108 
33.5 

.109 

34 

.201 

48.1 

.202 

48.3 

.203 

49.1 

.204 

49.1 

.205 

48.1 

.206 

49.9 

.207 

TAP 

1 

I . B83 

1.047 

1.038 

1.053 

1,059 

CP, 

1.078 

n 

1 ,073 

1.869 

1.873 

1.022 

1.822 

1.015 

1.023 

1.829 

1.029 

1.026 

a 

0 , 99t 

1.041 

0.984 

1.043 

1 . 01 * 

1.033 

1.063 

1.041 

0.918 

8.948 

1.821 

0.994 

1.011 

0.977 

0 . 99a 

1.025 

3 

1.030 

1.002 

1.031 

1.044 

1.039 

1.032 

1.081 

1.862 

1.067 

8.995 

1.822 

0.999 

1,014 

1.007 

1.007 

1.028 

a 

1.031 

0,994 

1.039 

l . 0«2 

1.043 

1,038 

1,057 

1.081 

1,075 

0,991 

1,013 

1.013 

1.012 

1 , 00b 

1.017 

1.022 

5 

0,995 

1.028 

1.050 

1.041 

1.032 

1 . 08 * 

1.069 

1.064 

1,079 

0.993 

1.015 

1,807 

1.025 

1.013 

1.018 

1,028 

6 

1.021 

1.043 

1.323 

1.333 

i . 03a 

1,059 

1.066 

1,064 

1.069 

1.001 

1.015 

1 . 00b 

1.032 

1 . 01b 

1.018 

1.027 

7 

1.018 

1 . 02b 

1.007 

1.352 

i , 03a 

1.053 

1.062 

1.062 

1.068 

1.012 

1,024 

1.007 

1.030 

1.020 

1.822 

1.019 

8 

1.017 

1.059 

1.325 

1 . 04b 

1.051 

1 . 05b 

1 . 06b 

1.063 

1.066 

1.023 

1.018 

1.009 

1.028 

1.012 

1.021 

1,030 

9 

1,017 

1 . 07b 

1.011 

1 . 0«4 

1.333 

1.063 

1.054 

1.065 

1.067 

1 . B26 

1.818 

1.014 

1.023 

1.021 

1.021 

1 .033 

10 

1.015 

1,048 

1.030 

1.043 

1.057 

1.040 

1.033 

1.065 

1.067 

1,814 

1.819 

1.018 

1.029 

1.015 

1.023 

1.038 

1 1 

1 . 00S 

1,014 

1.008 

1.041 

1.048 

1.027 

1.062 

1.065 

1.059 

1,026 

1.031 

1.813 

1.029 

1,024 

1.811 

l . a3 | 

>2 

0 . 99b 

1.007 

1.000 

1 . 0*1 

1.037 

1.032 

1.034 

1.062 

1.061 

1.023 

1.823 

1.818 

1.031 

1.022 

1.016 

1.034 

13 

0 . 99b 

1.037 

0.993 

1.016 

1.024 

1.032 

1,030 

1.059 

i . 85a 

1.012 

1.022 

1,813 

1,028 

1.017 

1.023 

1.058 

14 

0.993 

1.041 

0 . 99a 

1.014 

1.032 

1.029 

1.031 

1.853 

1,045 

1.014 

1.827 

1,022 

1.033 

1.019 

1.017 

1.056 

15 

0 . 99b 

1.022 

0,978 

1.032 

1.030 

1.030 

1.028 

1.063 

1.045 

1.019 

1.825 

1.020 

1.039 

1.825 

1.021 

1.040 

lb 

1.000 

1.058 

0.995 

1.028 

1.033 

l.a29 

1.047 

1.041 

1,043 

1.024 

1.027 

1.024 

1.035 

1 . 02b 

1.821 

1.037 

17 

0 , 99b 

1.080 

0.989 

1.031 

1.035 

1.030 

i . 0«4 

1 . 05b 

1.044 

1,827 

1.031 

1.023 

1 , 03b 

1.027 

1 . 02b 

1.053 

18 

1,012 

1.041 

0 . 99J 

1.027 

1.034 

1.017 

1.056 

1.038 

1.041 

1.018 

1.824 

1.823 

1 . 0*1 

1.026 

1.024 

1 . 04b 

19 

0.982 

1.040 

0 . 99a 

1.315 

1.039 

1 , 0*7 

1.048 

1 . 0*0 

1.047 

1 . 02b 

1.025 

1.887 

1 . 8«5 

1.032 

1.024 

1.047 

20 

1.210 

1 ,057 

1.005 

1.034 

1.038 

1,019 

1,053 

1.055 

1.037 

1 . 03b 

1.018 

1.019 

1.642 

1.041 

1.030 

1.851 

1 a) f 

21 

1.019 

1.037 

0.983 

1.009 

1.041 

1.020 

1.817 

1 , 0*2 

0,991 

1.032 

1.883 

1.021 

1.841 

1.021 

1.029 

1 • BC » 
a OhL 

22 

1.012 

0.881 

0.971 

0.905 

1,005 

1.024 

8.899 

1.040 

0,858 

1.027 

8.882 

1.021 

8.964 

0 . 99 | 

1.031 

* a 3 7 

23 

0,993 

1.027 

0.979 

1.018 

1.021 

1,037 

1.005 

1.046 

1.027 

1.004 

1,013 

1.013 

1.823 

1.028 

1.018 


24 

• 0.054 

• 0.032 

• 0 . 07a 

• 0.049 

• 0.032 

• 0.053 

• 0,064 

0.005 

• 0,045 

• 0.012 

• 8.812 

. 0.089 

• 0.004 

■ 0 . 80S 

• 0.009 

* 0.006 

25 

• 0.057 

• 0 . 05s 

• 0.074 

• 0,065 

. 0.049 

. 0,069 

• 0.081 

• 0,058 

• 0,068 

• 8.889 

• 8.821 

• 0,009 

• 0.014 

• 0,010 

• 0.023 

• 0801 ? 

2b 

0 . 9ba 

1,019 

0 . 9bb 

0 , 99b 

1 . 02b 

1,017 

1.003 

t.037 

1.016 

1.003 

8.990 

1.001 

1.818 

1.008 

1.089 

1 . 01 T 

27 

• 0.053 

• 0 . 03jS 

- 0 . 0/9 

• 0.044 

• 0.041 

• 0.060 

• 0,074 

• 0.847 

■ 0.067 

• 0,809 

• 8.882 

• 0.012 

8 . 80b 

• 0.010 

• 0.809 

« 0 * 00 ® 

28 

1.011 

1.044 

0.991 

1.021 

1.038 

1.055 

1.037 

1.868 

1.053 

1.838 

1 . 02 * 

1 . 03b 

1 , 0*4 

1 . 04 * 

1 . 0*4 

1 .045 

48 

• 0.057 

■ 0,045 

• 0.080 

• 0.057 

■ 0,039 

• 0,068 

• 0 , 07b 

• 8.843 

■ B , 06b 

• 8.808 

• 8.810 

• 0.007 

• 8.082 

• 0.005 

• 0 . 00 * 

* 0 .006 


NOTE: Three-digit descriptor .abc in file identifications of form TUNCAL.abc 

Designate: rake position by last digit (c ) , tunnel speed by first digit (a) [1,2, 3,4,5 nominally 35,50,70,100,140 mphl 
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APPENDIX XIX, Part A (contini 
Uncorrected calibration Data for TUNCAL fa b C 


FILE.abc 
Vel. ,mph 
TAP 

.208 
48. 6 

.209 

48.9 

.301 

70.3 

.302 

69.8 

.303 

70.2 

• 304 
69.4 

.305 

69.7 

.306 

69.1 

i 

1.823 

1.023 

1.008 

1.005 

1.007 

Cp 

1,010 

n 

1.813 

1 .009 

2 

0.986 

8 . 9ia 

0,947 

1.004 

0.982 

1.011 

0.976 

0.970 

3 

1.000 

1.025 

8 , 98 | 

1.804 

0.989 

1.013 

0.992 

0.985 

4 

1.089 

1.033 

0.973 

1.008 

0.992 

1.017 

0.997 

0.988 

s 

1.013 

1.024 

8.989 

1.082 

0.990 

1.020 

1.000 ‘ 

0.994 

6 

1.013 

1.037 

0.988 

1.015 

8.988 

1.018 

0.999 

0.996 

7 

1 ,015 

1.032 

0.991 

1.011 

0.996 

1.019 

1.802 

1.000 

8 

1 ,015 

1.033 

0.998 

1.010 

8.997 

1 . 01S 

1.001 

1.000 

9 

1 , 01b 

1.037 

0.999 

l.ais 

0.995 

1.010 

1.002 

1.081 

18 

1.017 

1.036 

1.082 

1.815 

1.003 

1.012 

1 . 085 

1.083 

11 

1,027 

1.830 

1.004 

1.018 

1.084 

1.016 

1.086 

1,084 

12 

1.026 

1.030 

1.304 

1.012 

1.086 

1,011 

1.085 

1.006 

13 

1,023 

1.035 

1.005 

1.023 

1.004 

1.016 

1.005 

1,009 

14 

1.023 

1.032 

1.012 

1.024 

1.009 

1.022 

1.006 

1.012 

IS 

1.016 

1.028 

1.088 

1.023 

1.009 

1.018 

1.013 

1.014 

lb 

1 ,03 5 

1.029 

1.017 

1.017 

1.000 

1.020 

1.811 

1.021 

IT 

1.027 

1.029 

1.012 

1.024 

1,084 

1,026 

1.020 

1.020 

18 

1.029 

1.024 

1.013 

1 ,025 

1.007 

1.026 

1,014 

1,019 

19 

1.035 

1,016 

1.013 

1.020 

1.007 

1.032 

1.019 

1.019 

28 

1.035 

1 .008 

1.014 

1.021 

1.006 

1,034 

1,018 

1.022 

21 

1.042 

1.012 

1.813 

1.013 

1.013 

1.024 

1.017 

1.021 

22 

1,026 

0.864 

1.018 

0,889 

1.012 

0,942 

0,989 

1.027 

23 

1.026 

1.024 

1.006 

1.081 

1.003 

1.014 

1.012 

1.811 

24 

0,040 

0.006 

• 0.004 

• 0.083 

• 0.084 

0,001 

• 0,081 

- 0.082 

f»5 

• 0,021 

• 0.018 

• 0.006 

• 0.014 

• 0,004 

• 0,014 

> 0,014 

• 0.017 

26 

1 .019 

1.014 

0.992 

0.992 

0.987 

0 . 99J 

0,990 

0,994 

27 

• 0.008 

• 0.813 

• 0.003 

0.083 

• 0.007 

0,009 

• 0,005 

• 0,003 

28 

1 .061 

1.040 

1.029 

1,022 

1.026 

| 1 . 0 J2 

1,036 

1.031 

48 

• 8.003 

• 0.007 

• 0.001 

< 0.080 

• 0.001 

1 0,080 

• 0.000 

0,000 


I ) 1 . . J '-.ZJ 11 --I)' — 1 — y — 1 
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APPENDIX III, Part A (continued) 

Uncorrected Calibration Data for TUNCAL.abc Piles, First (Model Station) Plane 


t 

FILE. ah c 
Vel. ,mph 

.405 

96.9 

.406 

99.5 

.407 

99.7 

.408 

100. 

.409 
100. . 

.501 .502 

125.8 124.5 

.503 

125. 

.504 

124.9 

.505 

124.3 

.506 

125. 

.507 

124.8 

.508 

125.1 

.509 

124.7 

TAP 

1 

2 

1 1 087 

1.003 

1.004 

1.082 

1.004 

cp n 

1.007 1.004 

1.009 

1.009 

1.010 

1.005 

1.003 

1.002 

1.004 

0.972 

8.958 

1.000 

0.978 

0,9lt 

0.920 

0.997 

0.957 

0.991 

0.963 

0.947 

0.987 

0,969 

0,908 

3 

0. 947 

0,980 

1.001 

0.993 

0,996 

0.966 

0.993 

0.977 

0,996 

0,982 

0.976 

0.994 

0.98J 

0.985 

4 

0.988 

0.986 

1,002 

0,993 

0.999 

0.966 

0.994 

0.982 

0.997 

0.979 

0,979 

0.993 

0,986 

0,986 

5 

8.992 

0,988 

1,002 

0,993 

1,003 

0.971 

0.994 

0,984 

0,999 

0,982 

0.981 

0.995 

0.986 

0.991 

& 

0.992 

0,988 

1,008 

0,995 

1.204 

0.979 

0.993 

0,984 

0,998 

0,984 

0.983 

0.998 

0.990 

0.994 

7 

0.99J 

8.991 

1.009 

0.995 

1.003 

0.979 

0.994 

0,987 

0.994 

0.980 

0.988 

0.998 

0.99J 

0.992 

a 

0.993 

0,994 

1.01S 

0.999 

1,004 

0.987 

0.993 

0.985 

0,998 

0.981 

0.988 

0.999 

0.993 

0.991 

9 

0.996 

0.996 

1,015 

0,999 

l,00« 

0.987 

0.991 

0.992 

0,998 

0.983 

0.991 

0,999 

0.996 

0.993 

ia 

0,999 

0.995 

1.015 

1.001 

1,007 

0.985 

0.992 

0.995 

0.997 

0.982 

0.992 

1.001 

1.802 

0.992 

XI 

12 

0,997 

0.995 

1,019 

1.002 

1,006 

0.991 

0.997 

0.997 

0.997 

0.990 

0.994 

1.001 

1,002 

0,995 

0,999 

0.999 

1.020 

1.005 

1.010 

0,990 

0.993 

0.998 

1.000 

0.991 

0.997 

1.001 

1.002 

0.99a 

13 

1.001 

1,004 

1.018 

1.007 

1.012 

0.994 

0,996 

1.000 

1.001 

0.992 

0,999 

1.002 

1.003 

0.996 

14 

1.003 

1,006 

1.021 

1.001 

1,009 

0.996 

0.996 

0.998 

1.002 

0.990 

0,998 

1.009 

1.004 

0.996 

15 

lb 

17 

1,004 

1.009 

1.023 

1,288 

1,012 

0.997 

0.994 

1.000 

0.999 

0,991 

0.999 

1.002 

1,003 

1.000 

1,004 

1.012 

1.026 

1.005 

1.011 

0.998 

0.997 

0,997 

1.001 

0.992 

1.001 

1.005 

1.001 

1.002 

1.008 

1.013 

1.029 

1,008 

1,007 

0.999 

1.000 

1.002 

1.084 

0,995 

1.001 

1.002 

1,000 

1.003 

18 

1,009 

1.014 

1.026 

1.004 

1.003 

1.000 

1.003 

1.000 

1.006 

0.994 

0.998 

1.003 

0.998 

1.004 

19 

2a 

21 

22 

23 

1,011 

1.013 

1,032 

1.007 

1,005 

1.000 

1.000 

1.001 

1.010 

0.996 

1.000 

1.005 

0.998 

1,005 

1,008 

1.019 

1.029 

1.007 

0.995 

0.999 

1.002 

1.001 

1,014 

1,000 

1.001 

1.003 

0,998 

0.997 

1,006 

1.018 

1.020 

1.004 

0,991 

1.001 

0,986 

1.001 

1.009 

0.999 

1.000 

0.999 

0,990 

0.986 

0.954 

1,022 

0.962 

0.986 

0.897 

1.001 

0.853 

1.002 

0,922 

1.001 

1.001 

0.919 

0,956 

0.894 

1.011 

1.005 

1.004 

1.005 

1.009 

1.009 

1.002 

1.012 

1.015 

1.010 

1.006 

1.005 

1.008 

1.007 

■ 24 

0.004 

0.001 

0,003 

0.044 

0.017 

0.005 

0.000 

0.006 

0.008 

0.005 

0,004 

0.006 

0.046 

0.020 

25 

26 

27 

28 

• 0.016 

•0.017 

•0.01 1 

•0.020 

•0,013 

•0.005 *0.012 

• 0,002 

• 0,013 

•0.015 

•0.015 

•0.014 

-0.017 

•0.012 

0.969 

0.982 

0.988 

0.982 

0.980 

0.980 

0.972 

0.977 

0.981 

0.974 

0.973 

0.97T 

0.974 

0,973 

•0.000 

0.003 

0.003 

0.002 

0.001 

0.009 

0.016 

0.005 

0.018 

0,003 

0.005 

0,005 

0.003 

0.004 

1.028 

1.003 

1.026 

1.019 

1.022 

1.011 

t.097 

1.013 

1.011 

1,009 

1.006 

1.008 

1.006 

1.008 

48 

0.001 

0.002 

0.002 

0.001 

0.081 

0.002 

0.002 

0.003 

0,003 

0.001 

0.002 

0.002 

0,002 

0.001 




Appendix III, Part B !• j <[> 

Pressure coefficient distribution across, tunnel section. Rake positions 1 and 2. 
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Appendix III, Part B (continued) 

Pressure coefficient distribution across tunnel section. Rake positions 3 and 4 0 

i,t.) i i i _j in ;_j :_J ' _) — > ~J 
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Appendix III, Part B (continued) 

Pressure coefficient distribution across tunnel section. Rake position 5. 



Pressure Parameter, Cp. 
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Appendix III, Part B (continued) 

Pressure coefficient distribution across tunnel section 
Rake position 8„ 



Pressure parameter, Cp 
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^ Appendix III, Part B (concluded) 

Pressure coefficient distribution across tunnel section. 
Rake position 9. 


J ■_ ) 1 J 1 — } ) ZiJ i • ■ i — 



Test Room Side 




Appendix III, Part C. Contours of Constant Pressure Coefficient Cp 
(continued) TUNCAL Phase of Test, Facing Upstream 

Wind speed 50 mph 
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Appendix III, Part C. Contours of Constant Pressure Coefficient Cp 
(continued) TUNCAL Phase of Test, Facing Upstream 

Wind speed 70 mph 
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Appendix III, Part C. Contours of Constant Pressure Coefficient Cp " 
(continued) TUNCAL Phase of Test, Facing Upstream n 

Wind speed 100 mph 

1 .... / 1 .. J . ) ) ) 


J .) 



Pressure Coefficient Cp 
Facing Upstream 



FILE. abc 
Vel . ,mph 
TAP 

1 

.241 

49.2 

.242 

49.4 

APPENDIX IV, Part A 

Uncorrected Calibration Data for TUNCA2.abc Files, Second (Downstream) 

.243 .245 .246 .247 ' .248 .441 ' .442 .443 .445 

49.7 48.7 49.7 49.5 49.2 100.1 98.8 100.1 100.1 

Plane 

.446 

100.1 

.447 

99.9 

1,009 

1,005 

1,011 

1.015 

0.995 

n 

1,009 

1.005 

1.008 

1,008 

0,999 

1,001 

1.003 

r ,003 

2 

0,885 

1.304 

0,939 

0,931 

0.901 

1,017 

0,932 

0.854 

0,998 

0.904 

0,949 

0.854 

0.992 

3 

0,949 

1.306 

0.998 

0,990 

0,960 

1.003 

0.973 

0.938 

0.997 

0,967 

0 , 9 b 7 

0.969 

0.994 

4 

0 . 9 b 6 

0,999 

1.004 

0.994 

0.975 

1,016 

0,987 

0.960 

1,004 

0,971 

0.975 

0.976 

0,998 

5 

0,987 

1.001 

1.003 

0,985 

0.979 

1,007 

0,990 

0,974 

1,005 

0,980 

0.978 

0,985 

1.000 

6 

0,983 

1,007 

1.007 

0,990 

0.984 

1,313 

0,995 

0,982 

1.304 

0,978 

0,976 

0,982 

1,006 

7 

1,002 

1.004 

1.302 

0,993 

0.986 

1,009 

0.995 

0.989 

1.013 

0.978 

0,970 

0,987 

1,001 

8 

0.988 

1,300 

1,306 

0.989 

0.982 

1.310 

0.989 

0.992 

1.310 

0.980 

0.975 

0.987 

1,010 

9 

0.995 

1.012 

1,310 

0,992 

0.986 

1.007 

0,992 

0,996 

1.011 

0,985 

0.978 

0.991 

1,039 

10 

0.557 

1.007 

1.332 

0.993 

0.986 

1,000 

0.987 

0.998 

1,308 

0,989 

0.987 

0,992 

1,009 

1 1 

0,994 

1.211 

1,312 

0,994 

0 . 9 B 5 

1.000 

0,993 

0,995 

1,009 

0.992 

0.987 

0.994 

1.313 

12 ’, 

0,987 

1.001 

1.312 

0.991 

0.981 

1.006 

' 0.985 

1,002 

1,012 

0.996 

0,989 

0,994 

1,015 


0,992 

1,301 

1.308 

0.995 

0.981 

1,002 

0.991 

1,003 

1,016 

0,995 

0,989 

1,000 

1 , 01 b 

14 

0,967 

1.303 

1,017 

0,985 

0.992 

0,993 

1.003 

1.000 

1.011 

0,994 

0,990 

1,300 

1.015 

15 

0,999 

0,990 

1,315 

0.992 

1,006 

0,993 

1.004 

1.003 

1,015 

0.994 

0,984 

1.305 

1.316 

lb 

0,997 

1,002 

1.308 

0,986 

0.66 2 

1,002 

0.998 

1.306 

1.021 

0.995 

0.991 

1.009 

1,018 

1 7 

1,000 

0.995 

1.210 

0.997 

0,629 

1.004 

0,993 

1,008 

1,020 

0.993 

0,987 

0.920 

1.320 : 

18 

1.301 

0,997 

1.013 

0,999 

0.717 

1.003 

0.998 

1.309 

1,020 

0,996 

0.996 

0,918 

1.321 

19 

1,201 

0,999 

1.310 

1.007 

0.792 

1.000 

0.996 

1.009 

1,017 

0,995 

0,996 

0.953 

1,320 j 

20 

1,002 

0.999 

1,009 

1,008 

0.853 

1.031 

1,001 

1.310 

1,008 

0.993 

0,997 

0.971 

1,317 

21 

0,994 

0.972 

1,206 

0.979 

0,899 

0,980 

0,986 

1,313 

0,993 

0,997 

0,977 

0.991 

1,016 . 

22 

1.000 

0,933 

1,006 

0,905 

0.928 

0.907 

0,991 

1,315 

0,893 

0,996 

0.924 

0.998 

0,959 

23 

• 0 , 3 i 2 

- 0,049 

- 0,330 

• 0,043 

" 0.051 

" 0,047 

0.012 

• 0.021 

• 0,025 

• 0.016 

• 0.018 

• 0.018 

• 0,020 

24 

0,994 

0,979 

0.995 

0.989 

0 . 9 b 2 

0,988 

0.994 

1.008 

1.000 

0,992 

0,998 

0.997 

1,007 

25 

0,977 

0,956 

0.972 

0.966 

0.931 

0,959 

0.960 

0.974 

0.969 

0,962 

0,964 

0,960 

0.971 

26 

* 0.032 

• 0.058 

• 0.044 

* 0,064 

• 0,065 

• 0.060 

- 0,067 

• 0.025 " 0.037 

• 0,034 

* 0.039 

* 0.039 

• 0.038 

27 

1,010 

0.991 

1.011 

1,004 

0.988 

1.001 

1.020 

1.320 

1,016 

1,009 

1,012 

1.014 

1,019 

28 

- 0.018 

• 0.026 

* 0,021 

" 0,035 

• 0.042 

• 0.035 * 0,336 

- 0,007 

- 0,003 

• 0,009 

• 0,007 " 0,009 

• 0.007 

48 

* 0.310 

- 0,027 

- 0,015 

• 0.029 

• 0,036 

- 0,027 

" 0,029 

• 0,003 * 0.003 

• 0,004 

• 0,002 " 0,003 

• 0,002 


NOTE: Three-digit descriptor .abc in file identifications of form TUNCA2.abc 

Designate: rake position by last digit (c) , tunnel speed by first digit (a) 

[1,2, 3, 4, 5 nominally 35,50,70,100,140 mph] 
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APPENDIX IV, Part A (continued) 

Uncorrected Calibration Data for TUNCA2.abc Files, Second (Downstream) Plane 


FILE.abc 
Vel. ,mph 
TAP 

i 

.448 

99.6 

.541 

124.6 

.542 

123.4 

.543 .545 

124.5 125.7 

.546 

124. 

.547 

122.8 

.548 

123.7 

1.002 

1.012 

1.009 

0,998 

1,000 

1,006 

1.005 

1.003 

2 

0 . 9 b 3 

0,824 

0,987 

0,904 

0.929 

0,883 

0,988 

0.960 

3 

0.990 

0.949 

0.991 

0,942 

0,954 

0,957 

0.992 

0.981 

a 

0.993 

0 , 9 b 0 

0.992 

3,960 

0,970 

0,974 

0,995 

0.985 

5 

0.994 

3 . 9 b 1 

0.992 

0 , 9 bb 

0.968 

0.977 

0.997 

0,987 

b 

0,993 

0 . 97 b 

0,994 

0 , 9 bb 

0.966 

0,980 

0,999 

0,907 

7 

0,994 

0 . 97 b 

0,999 

0.970 

0 , 9 b 9 

0,979 

1.000 

0,989 

8 

0 . 99 b 

0.983 

0,997 

0.969 

0,968 

0,985 

1.000 

0,992 

9 

0.997 

0.987 

0.996 

0.974 

0,974 

0,987 

1,004 

0.992 

12 

0.997 

0.988 

0.997 

0,979 

0,976 

0.987 

1.002 

0,998 

11 

0.998 

0.988 

0.997 

0.985 

0,975 

0.993 

1,000 

0,997 

12 

0.949 

0.993 

0,999 

0,982 

0,977 

0,994 

1.003 

1.002 

13 

1,002 

0,992 

0.999 

0,986 

0.979 

0,997 

1,000 

1,000 

14 

1,002 

0 , 99 b 

1,000 

0,987 

0.982 

0,999 

1,004 

1,200 

15 

1,003 

1.000 

1.001 

0,983 

0.977 

1,001 

0.999 

0,996 

lb 

1,005 

1.001 

1,004 

0,981 

0.978 

1,200 

0,999 

0.999 

17 

1.001 

1,000 

1.004 

0,978 

0,973 

1.003 

1.002 

0.997 

18 

0.999 

1.000 

1,006 

0,980 

0,980 

1,202 

1.001 

0,995 

19 

0,998 

0.999 

1.004 

0.981 

0,986 

1,001 

1,005 

0.995 

20 

0.997 

1.001 

0,997 

0,984 

0,990 

1,202 

1 , 003 

0.983 

21 

0,993 

1.001 

0,988 

0.986 

0,987 

1,203 

0,987 

0.975 

22 

0 , 9 b ) 

1.003 

0.936 

0,987 

0,948 

1.201 

0,959 

0.913 

23 

0,039 

- 0.012 

• 0,016 

• 0,008 

• 0,008 

• 0.009 

• 0,013 

0,044 

24 

0,998 

1.009 

1,004 

0.993 

0,998 

1.001 

1,004 

1.003 

25 

0.972 

0.959 

0.953 

0.948 

0,953 

0,960 

0.962 

0.959 

2 b 

• 0,041 

• 0.026 

• 0,037 

• 0.032 

- 0,041 

• 0,040 

• 0,036 

* 0,040 

27 

1,012 

1.008 

1,005 

0.999 

0,997 

1,000 

1.003 

1.200 

28 

• 0.009 

• 0.003 

0,001 

• 0,004 * 0,004 

• 0,006 

• 0.004 

* 0.006 

48 

• 0.003 

• 0.002 

• 0,002 

• 0,002 

• 0,003 

• 0.222 

• 0.001 

• 0.002 


J - 1 1 } 1 . ] I 1 i ... ) ) ) ! J : ) . ! J J J 



Pressure parameter, Cp. 
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Appendix IV, Part B 

Pressure coefficient distribution across tunnel section. 


Rake positions 1 and 2 



Pressure parameter, Cp 


1.05 


TUNCA2.243 


7 


50 mph 


.244 (estimate) 


100 mph 


,444 (estimate) 


125 mph 


.544 (estimate) 


Appendix IV, Part B (continued) 

Pressure coefficient distribution across tunnel section. Rake positions 3 and 4 (estimated) 



Pressure parameter, Cp 
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rt rr 
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50 mph 


1.05 
■ c 

1.00 
0.95 




• 445 100 mph 
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Appendix IV, Part B (continued) 

Pressure coefficient distribution across tunnel section. Rake position 5. 



Pressure parameter, Cp 







Appendix IV, Part B (continued) 

Pressure coefficient distribution across tunnel section. Rake positions 6 and 7. 




Pressure parameter, Cp 
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Appendix IV, Part B (concluded) 

Pressure coefficient distribution across tunnel section. Rake position 8. 
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APPENDIX B 


Adaptation of Two-Dimensional 
Transonic Analysis Code , TSFOIL , for 
Use in Modeling Adaptive Wall 
Test Sections 

May 23, 1983 

By: Jonathan Webb 
Advisor: Prof. E.M. Murman 



Introduction 


A difficulty encountered when performing wind tunnel tests 
is determining how to correct for interference caused by the 
walls of the test section. While all wind tunnel testing en- 
counters this difficulty, the effect of wall interference becomes 
most pronounced as the flow Mach number approaches unity. 

Shock waves forming on the aerodynamic surfaces propagate 
laterally instead of being swept downstream as in supersonic 
flow. This can choke the tunnel as in a convergent-divergent 
nozzle leading to gross errors in pressure as well as balance 
data. Numerous techniques have been developed over the years 
in an effort to combat this and other difficulties encountered 
with wind tunnel testing near Mach 1. Slotted or porous walls, 
correction parameters and enlarged test sections are among 
the most succesful techniques but they are all fixes to a prob- 
lem and are often dependent on engineering judgement when they 
are applied to an actual test. 

A technique currently being develeloped by NASA at Langley 
Research Center is an adaptation of a technique originally 
developed in England in the 19^0* s. For a two dimensional 
test, if the top and the bottom walls of the test section can 
be deformed so they are the same shape as a free air stream- 
line, all the effects of wall interference would be eliminated. 

There are a number of difficulties with this approach. 

The principle difficulty is, once the proper streamline shape 
has been found, the wall will never correspond exactly to the 
shape of the streamline. It is hoped the elastic deformations 



will be small enough so the wall held to a particular shape 
at only a few points will still reduce the interference 
considerably. 

The purpose of this project is to simulate numerically 
the effects of the elastic deformation of the wind tunnel wall 
held to a streamline-like shape. The wall itself is simulated 
by a finite element model and is allowed to deform under the 
pressure loading developing in the wind tunnel with an air- 
foil model present. A modified version of the transonic 
analysis program, TSFOIL, is then used to determine the resulting 
flow filed with the effects of the deformed wall included. 

The shape of the wall is updated during the course of the 
iterative solution of the flow field so the resulting solution 
correctly models the interaction between the shape of the wall 
and the flow field. Once a result from a particular operating 
condition has been obtained, the pressure distribution on the 
airfoil in the wind tunnel model can be compared with solutions 
generated with TSFOIL in free air mode. Any differences 
between the two solutions are primarily due to the interference 


effects of the deformed wall. 



The? Finite Element Model 


The? finite element model used to represent the wall is 
based on the principle of stationary, total potential energy 
(PSTPE) and the Rayliegh-Ritz procedure. Each segment of 
the wall between the jacks is treated as an individual element. 
A quintic, polynomial, bending only shape function is used to 
describe the deformation of the element. The von Karman plate 
formulation for the strain energy (see Reference 1) is used 
and is simplified by assuming the bending of the plate to be 
two dimensional. Use of the Rayliegh-Ritz procedure and the 
PSTPE result in a system of linear, algebraic equations which 
are solved by point relaxation. The resulting shape of the 
wall is then included in the flow field calculation. 

For two dimensional, plate bending with no axial loading, 
the expression for the internal strain energy is: 

u fe = h * j D(w") 2 dx (1) 


where: 

u fe is the bending strain energy density per unit width 
w is the bending deflection 

' denotes differentiation with respect to the axial 
coordinate, x 


D = 


E tf 
12(l-v 2 ) 


E is Young's modulus = 2.8 x 10^ psi 
t is the plat thickness 





v is Poisson's ratio = 0.3 

By equating the strain energy to the work done on the 
element by the externally applied loads, the PSTPE can be 
applied and the unknown constants in the shape function can 
be found. 

The work done by the externally applied loads is given by: 

u e = “ f P z * w dx (2) 

where: 

u Q is the external work density per unit width 
6 

p .is the applied load per unit area in the bending 
direction 

The Rayliegh-Ritz procedure uses an assumed function 
which will describe the displacements of the element. For 
this model, a quintic polynomial shape function is used. 

Phis shgpe is: 

o o h, c 

w = a^ + a g x + a^x + a^x J + aye +• a^x^ (3) 

where the a^'s are unknown constants to be determined. 

It is desirable to describe the deformation fo the element 
in terms of displacement quantities which are defined at the 
ends of the element. These generalized displacements are 
described in Figure 1. 


• -.v. **F*«*6G-' •••• " v - ‘ -y - • ,,.u .. -•v, eOTiS r i .., 



The generalized displacements (q's) can be found in terms 
of the coefficients of the shape; function by using the following 
transf ormation . 


(a ) = J (q ) 
n' VH n 


where : 


( a^ ) and (q n ) are column matrices of q n and a r 

J = 1 00000 

0 1 0 0 0 0 

0 0 1 0 0 0 

-10 -6 -3 10 -4 1 

l 3 l 2 21 l 3 l 2 21 


li 3 J zkl 1 zk 

37 l 3 21 2 27 2 ? l 2 

zk zl zk k zl _1_ 

l 5 37 21 3 l 5 27 21 3 

If the shape function written in terms of the generalized 
displacements (using eq. 4) is introduced into the expression 
for the strain energy, the following expression results. 


u >, ■ i 


J D (a n ) T <x n )"(x n )" T (a n ) dx 


The integral can be evaluated to givei 






E* = 


0 

0 

0 

0 

0 0 





0 

0 

0 

0 

0 0 





0 

0 

41 

61 2 

81 3 101 4 





0 

0 

61 2 

121 3 

181 4 241 3 

(5) 
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81 3 

181 4 

l44l 3 401 6 





0 

0 
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4 241 5 
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401° 400r 






» 
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where 1 is 

the length 

of the segment 

• 



The transformation given by 

eq . 

4 and the element 

stiff 

ness matrix 

defined ast 







K = : 

D J T 

E* J 




(6) 


where : 









K = 

120 

6o 

_2_ 

-120 

60 

z2_ 



D 

737 

CM 
i — i 

71 

71 3 

7? 

71 




60 

1^2 

-24 

-60 

108 

-4 




71 2 

351 

35 

7i 2 

351 

35 




1_ 

-24 

Z61 

zl 

22- 

1 _ 




71 

35 

70 

71 

35 

70 




-120 

-60 

zl_ 
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-60 
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7 1 3 
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71 

71 3 

71 2 

71 




60 

108 

21 

-60 

1£2 

-11 
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351 

35 

71 2 

351 
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1_ 


-11 
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35 

70 

71 
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The generalized displacements can be introduced into the 


expression for the external work and used to find the gen- 
eralized loads. Defining the displacement, w, in terms of 
the generalized displacements yields* 

w = (D*) T (q n ) (7) 

where i 

(D*) = 1 - 10(x/l) 3 + 15(x/l) 4 - 6(x/l) 5 

x - 6l(x/l) 3 + 8l(x/l) 4 - 31(x/l) 5 
x 2 -31 2 (x/1) 3 + 1.51 2 (x/1) 4 - il 2 (x/l) 5 

10(x/l) 3 - 15(x/l) 4 + 6(x/l) 5 

-4l(x/l) 3 + 71(x/l) 4 - 31(x/l) 5 

|l 2 (x/l) 3 - l 2 (x/l) 4 + |l 2 (x/l) 5 

Introducing these expressions into the integral for the 
external work gives: 

u 9 = -Jp z (D*)' r (q n ) dx 

or 

u e = - (q n ) T (Q) (8) 

where 

(Q) = J p z (D* ) dx (9) 



Applying the PSTPE to the individual element yields: 

K (q n ) = (Q) (10) 

These expressions are only valid for an individual elemnt. 
To model the entire wall, an assemblage of these elements is 
required. By equating certain displacement quantities at i 
junctions between elements, a single expression of the same 
form as eq. 10, valid for the entire wall, can be determined. 
For the quintic shape function used, w' and w" are required 
to be continuous at the element junctions. This is the same 
as requiring: 


n_ n+1 

q 5 - q 2 


and 


n _ n+1 

q 6 " q 3 


where the superscript n is the number of the individual 
element. 

The generalized loads must be assembled in a somewhat 
different way. The load associated with qg» q^ . q^ or q^ 
must be added to the load from the adjoining element indent- 
ified by the equalities above. 

To complete the assemblage of the local elements, the 
boundary conditions must be included in the globil equations. 
At the upstream end of the wall, w’ is set to 0. At the 
downstream end of the wall, w" is set to 0. Since the 



element junctions are also the jack locations, w is prescribed 
at each node. Elements of the global stiffness matrix corre- 
sponding to the generalized displacements known from the 
boundary conditions are removed. Known elements effecting 
unknown displacements are included in the terms for the global 
generalized loads. This is so the only components of the 
global stiffness matrix remaining are those associated with 
the unknown, generalized desplacements . 

Once the global equation has been determined,, the resulting 
system of linear, algebraic equations is inverted using point, 
successive, over relaxation. Tis is not a method normally 
used for solving this type of system but an S.O.R. subroutine 
was readily available and could be used with very little 
modification. Reference 2 outlines the general procedure used. 

Once the generalized displacements have been determeined 
from the relaxation procedure, the deformation of the wall can 
be easily found using the J transformation matrix given by 
eq, 4. Wall displacements are determined at mesh points used 
during the current iteration in TSFOIL and the slopes of the 
wall at those mesh points is determined by using the cubic 
spline subroutine included in TSFOIL. The wall slopes are 
then used to modify the wall boundary condition and the flow 
field calculation proceeds. 



Modifications to TSFOIL 


The transonic, two dimensional, small disturbance analysis 
program, TSFOIL, written by Earll Murman, Frank Bailey and 
Margaret Johnson and documanted by Stephen Stahara in NASA 
GR-3064, was modified for use in this study. The modifications 
allowed the use of a variable geometry wall when specifying 
the boundaty conditions with the shape of the wall predeter- 
mined. The downstream boundary condition was also modified 
to better simulate the test section being modeled. 

In the original version of the program, a zero gradient 
condition was used at the downstream boundary of the com- 
putational domain. In the version used for this study, the 
downs tream condition was based on the one dimensional, isen- 
tropic, gas flow relation! 

du _ -dA * _l_p (11) 

u A l-»r 

where: 

u is the exit flow velocity 

A is the exit cross sectional area 

M is the exot Mach number 

Equation 11 can be differentiated with respect to the 
axial coordinate, x, and, upon introducing the perturbation 
velocity potential, can be written ass 


• SEKSMWf.,. 







( 12 ) 



-dA * 1 » 

dx A(l-M 2 ) 


This expression can be written in difference form for a 
non-uniform mesh as: 


n = * (x in - x i-i> (x i+i - x i) * B 


x i+i - x i 
x i " x i-l 


(13) 


(<v . - <fl. . ) + <$. 

V4 i Y i-l' i 


where B - 


dA * - 1 a 

dx A(l-M 2 ) 


Equation 13 is used as an extrapolation formula at the 
downstream edge fif the computational mesh. 

' The other modification to the boundary conditions involved' 
the treatment of the walls of the test section. The rriginal 
version of the program allowed only straight walled test 
sections to be modeled. For a wall with non-zero slope, 
the boundary condition must become: 


d£ = &L (14) 

dy dx 

where ^ is the slope of the wall. 

Since it is generally sufficient to include the boundary 
conditions to first order, eq. 14 can be written in finite 
differenc form as: 



" ^w+1 


(15) 


- y w } ^w+1 

where the w subscript is the index at the wall. This form is 
used for the upper wall. For the lower wall, the direction 
of indexing is reversed. 

Equation 15 is used as an extrapolation formula in the 
same manner as eq. 13 was used. 

It was also necessary to modify a number of the subroutines 
in TSF0I1 to interface the wall calculation with the flow field 
solution during the iteration process. The following is a list 
of the subroutines modified and a brief description of the 
modifications made to each. 

TSFOIL was modified to call the subroutines generating 
the global stiffness matrix for the wall deformation analysis. 

The subroutine which generates the wall shape during a free 
air calculation is also called by TSFOIL. 

READIN was modified to include the new input parameters 
in the new MIT/JCF VAX namelist statements. 

BLOCKDATA was modified to include all of the new common 
blocks as well as default values for the input variables. 

RESET was modified to include the boundary condition 
modifications described above. The vortex-doublet solution 
for the far field boundary was also removed from this subroutine. 
Since the computational mesh had to be extended to include 
the entire variable geometry wall, it was decided hhis 
portion of TSFOIL ' s original form would not be needed. 



Since this flow field condition; interacts with the program in 
many ways, it is not known how removing this section effects 
the remainder of the program. 

FARFLD was modified to treat the variable geometry wall 
far field condition the same as the solid wall case of the 
original program. 

SOLVE was modified to call the wall deformation routine. 

The follow, ng are new subroutines added for this study; 

STIFFMAT generates the global stiffness matrix for the 
strucural solution. It calls subroutine KLOCAL which calc- 
ulates the elements of the local element stiffness matrix 
needed for the global matrix. 

WALLV is the controlling program for the deformation 
routine. It determines the pressure distribution aiding the 
wall form the perturbation velocity distribution, it calls 
the deformation routine and it determines the correctly 
scaled wall slopes to be used in the wall boundary condition. 

DEFORM is the structural deformation analysis. This 
subroutine determines the global, generalized loads from the 
pressure distribution along the wall and the local wall 
geometry. Functions Q and BC determine the distribution 
factors for the pressure loads and the geometric boundary 
conditions respectively. The relaxation subroutine (sub- 
routine RELAX) is also called from this subroutine. When 
the generalized displacements have been determined, DEFORM 
also determines the actual wall shape so subroutine WALLV 



can determine the wall slopes as needed far the boundary- 
conditions. 

STREAM is the subroutine used to generate the streamline 
shape at the height of the tunnel walls. It is activated at 
the end of a free air case and can be used to generate the 
streamline shape internally for the wall deformation routine . 



Progress to Date 


In the original formulation of the deformation solution, 
a simplified form of the equilibrium plate equation was 
solved over segments of the wall in a manner similar to the 
finite element solution described previoulsy. This methos 
alos led to a system of linear, algebraic equations for a set 
of unknown coefficients in a polynomial describing the def- 
lection of the wall. This method was coded bjit a solution was 
never obtained from this procedure. It is not known whether 
the resulting system of equations was numerically ill-conditioned 
of if there were still undetected code errors. 

The finite element procedure was introduced when no 
further progress was being made with the equilibrium method. 

The finite element methos has not yielded a solution either; 
however, progress was still being made when work was stopped. 

The boundary conditions for the variable wall shape have 
been tested by using the streamline input to WALLV but bypassing 
the deformation analysis in DEFORM. The rleaxation solution 
seems to progress properly in TSFOIL but the results from 
the test case used have not been examined in detail nor have 
numerous test cases been run. 

All of the cases for the variable wall geometry encounter 
fatal run-time errors in subroutine PRTWAL. The version 
being used in this study had been modified for a previous 
project and the error seems to be connected with those mod- 
ifications. I do not believe using the original version of 
PRTWAL would alleviate the problem. 



New Input Parameters 


BCTYPE 

SLINEU 

SLINEL 


*** Note *** 


Integer. Set to 7 to activate the var- 
iable geometry wall solution. 

Real. Height of the wall from the centerline 
of the test section. Heights are positive 
for the upper wall and negative for the 
Shower wall. Input one value for each point 
in the input mesh in order from upstream to 
downstream, formalized to nominal 6" 
airfoil chord. Default i 1.1 and -1.1 
If a free air case is run with STREAMGEN = 

.TRUE. , when allowed to run to the fine mesh, 
subroutine STREAM will automatically set up this 
data if the variable geometry wall is run without 
terminating the program between cases. 


FREETEMP 

FREEPRES 


NWALL 


WALLEN D 


Free stream static temperature and pressure. 
Must be input in degrees Rankine and psf . 
Default: standard temperature and pressure. 
Integer. Number of flow field iteration 
performed before the wall shape is updated. 
Default: 10 

Logical. If .TRUE., the downstream end of 
the wall is set by the streamline input. 

If .FALSE., the last jack setting is set 
to the opening height given in the blue- 
prints fo r the wall. Default: .TRUE. 



STREAMGEN Logical. If .TRUE., a streamline is gen- 


WALLDEF 


erated starting at +H and -H as defined 
by input. Default* .FALSE. 

Logical. If .TRUE., the wall deformation 
analysis is called during a variable geo- 
metry wall case. Default: .TRUE. 
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Introduction and Summary 

This report describes a program for calculating a C-type mesh around airfoil like shapes. The 
Jameson/Caughey approach is used: a parabolic transformation coupled with a shearing 

transformation. The algebraic algorithm is capable of efficiently generating nearly orthogonal grids. 

A high degree of grid control is possible. The user may specify grid boundaries, number of grid 
lines, and location of (and spacing in) trailing edge and leading edge fine mesh areas. The 
capability of embedding fine mesh regions, for use with new adaptive grid techniques, is being 
developed. 

Grids generated by the program have been used in Euler flow calculations by W. Usab [5]. 
Results superior to results calculated on previous O-type grids were obtained. Specifically, 
calculations converged faster using C-grids rather than O-grids, total pressure loss spikes at the 
trailing edge of the airfoil were eliminated, and the Ni method converged with zero artificial 
smoothing for a subcritical case (resulting overall total pressure loss was then nearly zero). These 
improvements were attributed to higher grid orthogonality, especially at the trailing edge. 

This report is meant to serve as an outline of the C-mesh generator algorithm and its use. The 
program itself, is fairly straightfoward. Roughly half of the 800 code lines are comment lines. 


’"MB 
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1. Technical Background 

This discussion follows closely that found in [2]. Consider the physical space to be described in 
Cartesian coordinates, where x and y represent the streamwise and vertical coordinates respectively. 
We introduce a singular point just inside the leading edge of the airfoil. This singular point will 
be the locus of branch points in subsequent transformations that unwrap the airfoil to a shallow 
bump. The location is chosen to make the bump as smooth as possible. 1 We represent the 
singular line as: 



y =■• y s 

and define 

x = x - x s 

y = y - y s 

We then introduce the conformal mapping 
Xj + Y x = (x + y) 1 / 2 

which, after equating real and imaginary parts gives us two families of orthogonal parabolas, 
Y - Y = x 

and 

2X x * Yj = y. 

This transformation maps the entire airfoil to a shallow bump near Yj = 0 (see Fig. 1). 


If the nose of the airfoil can be represented by a parabola, the optimal branch point is the focal point. 
Otherwise, trying several values of x ? close to .005(chord) has quickly given an optimum for the geometries tested to 
date. 


swtfsr • ' 



Fig lb. Conformal Transformation with (x s ,y s ) = (0.006,0.). 

Figure 1. 

(Absolute values of Yj are used to allow more accurate spline fit.) 


Defining the height of this bump as 


Y, == S(X,) 


the final shearing transformation, 
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* - X l 

r, == Yj - S(X,) 

reduces the airfoil to t) = 0. Clearly, the orthogonality of the final grid depends on the height 
and smoothness of this bump which are functions, of the thickness distribution along the chord of 
the airfoil. 

1.1 Fine Mesh Regions 2 

Finally, we introduce a stretching transformation in 

S = Mj * y t 

between the trailing and leading edges, and 

* =: M 2 * y 2 

(Mj and M 2 are constants) from the trailing edge to the farfield. 3 y x and y 2 were specified as 
fifth and second order polynomials, 

=== a s^’ 5 + + - + a o 

and 

y 2 — l) 2^’ 2 + b i^.’ + b o’ 

respectively. We merely are stipulating that for any constant increase dy t or dy 2 we generate a 
proportional number of £. lines; therefor, the steeper the polynomial is, the higher is the number of 
£ lines. By varying the constants a Q - a 5 and b Q - b 2 we can control the spacing between 
subsequent £ lines. This is essential for designing grids with high resolution at the leading and 
trailing edges. 

The fifth order polynomial used between the trailing and leading edges, is uniquely determined 
by six boundary conditions. The start and end points (location of leading and trailing edge) 


~T|ie algorithm described in t his section was developed and programmed by W. llsab. 




Tjiis scheme has been applied only to f. lines here, but is easily extensible to n lines. 



determine two constants. The remaining four constants are determined by the slopes and 
curvatures at the end points. Thus, for example, if we specify a high slope at the leading edge 
and low slope at the trailing edge we will generate grid lines clustered in the nose area and sparse 
approaching the trailing edge. 

Three conditions must be specified for the parabola. The two endpoints, farfield and trailing 
edge locations, account for two of these. As the third condition, we require that the spacing 
between the two £ lines immediately after the T.E. equal the spacing of the two £ lines just before 
the T.E.. 

A sample grid, the grid used by W. Usab in the calculations mentioned in the introduction, is 
shown in Fig. 2, below. 



Fig. 2: 

17x79 Mesh around NACA 0012 airfoil with grid line packing in 
T.E. and L.E. areas. 
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1.2 Embedded Grids 

Embedded grids are another means of achieving high resolution near critical areas. The present 
program is capable of embedding any number of grids (M) at a given grid level (L) at arbitrary 
locations. We define a level as a collection of grids with equal spacing, Aq and A£, in 

computational space. 4 

The user supplies the chordwise and normal locations of the desired embedded region in physical 
space, as well as the spacing at each level, which must be an integral multiple of the global 
spacing. The program then calculates the locations in computational space that correspond to the 
physical chordwise and normal locations. 

The program first calculates the global mesh, saving the location of the grid nodes in 
computational space. Starting from the top level L, our embedding algorithm then steps through 
the global nodes in computational space in the positive £ direction, comparing the value of (; at 
each node to the user specified £ t t and £ d . After it locates the nodes that most nearly 

correspond to the user specified boundaries of the first mesh M x at L it simply divides that 

embedded region into (£ end - £ start )/A£ node points and labels each with the appropriate value of 
i;. It then does a similar search in the q direction. Having finished Mj at L, it repeats the 
procedure with the other meshes M 2 3 (if any) at L and then proceeds to lower and lower levels. 

Some sample plots are shown in Figs. 3 & 4. Presently, the routine is being refined and adapted 

for use with the pointer system of W. Usab. 


’Note, however, that A«) = A£ is not required. 
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Fig. 3: 

13x71 Grid with 2 levels of embbeded meshes at .7 chord around 
Korn & Garabedian airfoil. Level 1 has global/4 spacing (in both 
t) and £). Level 2 has global/2 spacing. 
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2. Program Capabilities 

The program generates a C-type mesh using a square root transformation coupled with a 
shearing transformation. For airfoil like shapes, the mesh is nearly orthogonal. 

The program prompts the user for: 

1. Grid outer radius (horizontal distance from mid-chord to front of grid). 

| 

-r 2. Grid extension (distance from trailing edge to end of grid). 

3. Grid spacing (in normal and tangential directions). 

;■ 

4. Grid line packing parameters (Note: This section is still being developed, so the 

operations haven’t been totally "mechanized” yet.) 

"j 5. Locations of up to 10 embedded 'meshes per level. (Any number of levels of embedded 

meshes are also possible, but the highest level meshes must be input first.) 

I 

The user may also input the singular point (x g , y s ). A plotting subroutine provides optional 
plots of the mesh. 

The program was written in Fortran IV and compiled and run on a VAX 750. Interactive 
graphics terminals were used at all stages of the grid design. 







-- JgfJ* J?' 
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3. Program Structure 

A flowchart of the computer program is given in Fig. 5, below. 

C-Mesh Generator Flowchart 


INPUT SECTION 

Read airfoil geometry i x,y ) 

Solve mapping equations to go to 
X t , Yj plane. 

Fit cubic spline to airfoil in 
Xj. Yj plane. 

Read boundaries of grid. 

Calculate corresponding values of £ ti 
with Newton routine. Read locations 
of fine and embedded grid areas and 
degree of fineness of grid. 

SOLUTION ALGORITHM 

Solve for £ and tj. 

Interpolate with the cubic spline 
function to find shearing term. 

Solve for Xj and Yj. 

Solve for x and y. 

Increment £ to £ . 

Increment q to tj bu . 

OUTPUT SECTION 

Write and plot results. 

Figure 5. 

The algorithm fits a cubic spline to the first transformed coordinates of the airfoil. It is used to 
interpolate to find values of Y — S(Xj), the shearing term. A Newton iteration routine is called 
whenever values of £ and ti are needed from stipulated x and y (when setting the limits of the 
grid, for ex.). 







4. Conclusions 


An algorithm for calculating C-type meshes around airfoils using the method of 

Jameson/Caughey is described. Generated grids are shown to .be very nearly orthogonal resulting 
in significant improvements in flow calculations compared to previous results. 

The program is very flexible. Fine mesh regions at trailing and leading edge regions may be 
defined conveniently. Any number of embedded grids with arbitrary spacing may be located 

within the global grid. Global grid characteristics such as outer radius, number of grid lines, etc., 
are modified easily. Several sample grids are given. 
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